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When polarized particles are accelerated in a synchrotron, the spin precession can be periodically 
driven by Fourier components of the electromagnetic fields through which the particles travel. This 
leads to resonant perturbations when the spin-precession frequency is close to a linear combination 
of the orbital frequencies. When such resonance conditions are crossed, partial depolarization or 
spin flip can occur. The amount of polarization that survives after resonance crossing is a function of 
the resonance strength and the crossing speed. This function is commonly called the Froissart-Stora 
formula. It is very useful for predicting the amount of polarization after an acceleration cycle of a 
synchrotron or for computing the required speed of the acceleration cycle to maintain a required 
amount of polarization. However, the resonance strength could in general only be computed for 
first-order resonances and for synchrotron sidebands. When Siberian Snakes adjust the spin tune 
to be 2, as is required for high energy accelerators, first-order resonances do not appear and higher- 
order resonances become dominant. Here we will introduce the strength of a higher-order spin-orbit 
resonance, and also present an efficient method of computing it. Several tracking examples will 
show that the so computed resonance strength can indeed be used in the Froissart-Stora formula. 
HERA-p is used for these examples which demonstrate that our results are very relevant for existing 
accelerators. 



I. INTRODUCTION 



In this paper we want to introduce the strength of 
higher-order spin orbit resonances which we want to use 
in the Froissart-Stora formula to compute how much po- 
larization is lost when a resonance is crossed. For first- 
order resonances the definition and computation of the 
resonance strength is relatively simple 0, , for higher- 
order resonances it is much more elaborate. We will need 
to use the invariant spin field, also called the n-axis [3|, 
the amplitude dependent spin tune and the periodic coor- 
dinate system over phase space that determines the spin 
tune 0. These concepts are therefore quickly reviewed 
in this introduction. 

While a polarized particle moves along the azimuth 
9 = of the storage ring's closed orbit with path length 
I and total length L, its semi-classical spin precesses ac- 
cording to the T-BMT equation 0, 



ae 



S = fi o (0) x S . 



(1) 



The spin direction that is periodic after one turn is re- 
ferred to as rio(9). If the spin has any other direction, 
it precesses around Hq . The numbers of precessions that 
occur during one turn is referred to as the closed orbit 
spin tune v . To describe the precession, a right handed 
system of orthogonal unit vectors (m, I, Ho) is introduced 
for any azimuth. The two vectors 1(6) and rh(9) precess 
around no according to the T-BMT equation so that they 
would have rotated v§ times after one turn. However 
a precession is added that continuously winds back i/q 
precessions. These vectors are therefore periodic in the 



azimuth and 



d9 



d9 



I = 



(fi - u n ) x fn , 
(fio - v n ) x / . 



(2) 
(3) 



Since particles on the closed orbit have spins that pre- 
cess around no, the product S3 = S(9) ■fto{9) is an invari- 
ant, i.e. does not depend on 9. It can be shown that it is 
also an adiabatic invariant i- e - ^ hardly changes 

when parameters of the system, like the storage energy, 
are slowly changed. 

This concept of an invariant spin direction, a spin tune, 
a periodic system of unit vectors and an adiabatic invari- 
ant can be extended to particles that do not move on 
the closed orbit but oscillate around this orbit and whose 
motion is thus described by phase space trajectories z(9). 
The T-BMT equation for spin motion then depends on 
the phase space trajectory 



d6 



s = n(z, (?) x s 



(4) 



If the vector field f(z,9) with \ f\ — 1 describes the 
spin distribution in a particle beam, it is called a spin 
field and satisfies the T-BMT equation 



d6 



f(g(6),6) = n&6)xf 



(5) 



A special spin field that is periodic from turn to turn is 
called the invariant spin field n, 



n{z, 9 + 2tt) = n(z, 9) 



(6) 
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Particles that travel along the trajectory z(9) have 
spins that precess around n{z{9),9). Describing this pre- 
cession and even the number of precessions in one turn 
starting at z(9q) is not trivial, since the particle has a 
new phase space point z(9q + 2n) after one turn. An 
orthogonal set of unit vectors (ui,U2,n) has to be de- 
fined for each phase space point and for each azimuth to 
determine spin precession angles. 

If the unit vectors u\ and u 2 would satisfy the T-BMT 
equation along each phase space trajectory starting at 
z*i and ending at Zf after one turn, these vectors would 
precess around n and after one turn Ui(zf, 9q + 2tt) would 
have some angle 2ttv{zq) with respect to the initial unit 
vectors Ui(zf,9o) at the same phase space point. 

The rotation angle 2-kv is not well defined, since the 
direction of the it, before and after the turn is only re- 
quired to be perpendicular to ft, but has a free angular 
orientation in the orthogonal plain. This free orientation 
for each phase s pac e point can (under certain general 
conditions HGUlni) be chose n to make the number of 
rotations v independent of the orbital phase variables 
It then only depends on the amplitudes J of the orbital 
motion and is therefore called the amplitude dependent 
spin tune v(J)- 

To obtain a periodic set of unit vectors, the described 
precession of the unit vectors is again augmented by con- 
tinuously winding back v spin precessions during one 
turn, 

^-Ui(z(9),9) = [Q(z, 9) - v(J)n(z, 9)] x u t . (7) 

Since spins precess around ft, the product 
Js = S(9) ■n(z(9),9) is an invariant of motion, i.e. it 
does not change with 9. It can be sho wn, however, that 
it is also an adiabatic invariant ELSGj, i.e. it hardly 
changes when system parameters like the storage energy 
change sufficiently slowly. This has strong implications. 
When a beam is polarized parallel to the invariant spin 
field n(z, Ei) at some initial energy Ei and the storage 
energy is increased slowly, the beam will be polarized 
parallel to n(z, Ef) at the final energy Ef. 

This is a very important property since a beam in such 
a polarization state will have the average polarization 
P\im —< n > after acceleration, which can be large even 
if this average polarization is small at intermediate ener- 
gies. 

II. THE SINGLE RESONANCE MODEL (SRM) 
A. Fourier Expansion of Spin Perturbations 

The quantities ft, v, u\, v,2 and Js will be computed 
for an analytically solvable model and the adiabatic in- 
variance will be illustrated by letting a parameter of this 
model change. Since this model leads to the Froissart- 
Stora formula, a comparison of its equations with the 



equations of general spin dynamics leads to the intro- 
duction of higher-order resonance strengths that can be 
used in the Froissart-Stora formula. 

The spin precession vector for particles which oscil- 
late around the closed orbit can be decomposed into the 
closed-orbit contribution VLq and a part Co due to the 
particles' oscillations, (l(z,9) — (lo(9) + Co(z,9). In the 
(m, I, rio) system we write 

S = Sim + S2I + S3W0 , Co = ujim + LU2I + • (8) 

With the complex notation s = si +1S2 and u> = io\ -\-\l)2<, 
the equation of spin motion is 

n x S = m ^ s i + l ^Q S2 + ~j^ s 3 + (^0 - ^o«o) x S (9) 

and the equation of motion for s is obtained by mul- 
tiplication with rh + il, and taking into account that 

■^-s = i(vQ + uj 3 )s - - \s\ 2 . (10) 

CLu 

In a coordinate system that rotates by vq9, this equation 
becomes 

So = e-' Wo6 s , ^§0 = ia^o-ie-^Wl-N 2 • (11) 

Spin precession on the closed orbit (ui = 0) leads to a 
constant §0 due to the left equation. The right equa- 
tion describes additional precessions due to phase space 
motion. 

If the motion in phase space can be transformed 
to action-angle variables, the spin precession vector 
uj(J, <&, 9) for particles which oscillate around the closed 
orbit is a 27r-periodic function of $ and 9. The Fourier 
spectrum of u>(J, &o + Q9, 9) has frequencies k = jo+j-Q 
where the jk are integers and Q describes the tunes of 
synchrotron and betatron oscillations. The integer con- 
tributions jo are due to the 2ir periodicity of Co in 9 and 
give rise to so-called imperfection resonances. The con- 
tributions j ■ Q of integer multiples of the orbit tunes are 
due to the 2ir periodicity of Co in the orbital phases $fe 
and give rise to so-called intrinsic resonances [j]. When 
one of the Fourier frequencies is nearly in resonance with 
Uo, one component of e~ lp ° e uj is nearly constant. Then 
it can be a good approximation to drop all other Fourier 
components since their influence on spin motion can av- 
erage to zero so that they are in effect less dominant. 
This is referred to as the single resonance approxima- 
tion. Note that this approximation can only be good 
when the domains of influence of individual resonances 
are well separated. This model corresponds to the ro- 
tating field approximation often used to discuss spin res- 
onance in solid state physics 01 ■ Note also that for a 
conventional flat ring, the first-order resonances due to 
vertical motion dominate and therefore the Fourier com- 
ponents with frequencies k = jo ± Q y are often of most 
interest. 
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The amplitude of a single Fourier contribution is some- 
times called the resonance strength. This is misleading 
since generally it cannot be used in the Froissart-Stora 
formula. The fact that the Fourier component is not 
the resonance strength manifest itself clearly in models 
where where uj is linear and has only first-order Fourier 
components, i.e. those with J\=i life I = 1- Such a uj can 
lead to depolarization or spin flip at first-order resonances 
but also at higher-order resonances 

EGHHEI. The 

strength of these resonances that might be use in the 
Froissart-Stora formula can clearly not be determined 
by the higher order Fourier coefficients, i.e. those where 
life I > 1j since those are zero. In fact all examples 
of higher-order resonances that will be shown in this pa- 
per were computed for such a linear model of HERA-p 
with Siberian Snakes |3Cj . 

A higher order resonance can thus be created either by 
a higher-order Fourier component or by feed-up of lower 
order components. Such a feed-up can occur due to the 
inherent non-commutativity of three dimensional rota- 
tions or equivalently due to the nonlincarity of the map- 
ping from the unit sphere to the complex plane which 
gives rise to the square root term in the equation of mo- 
tion HI U| l • Obtaining a resonance strength e K that can be 
used to describe depolarization therefore has to include 
all these feed-up effects. Before the following investi- 
gations it was not clear whether a Froissart-Stora for- 
mula with some resonance strength e K could be applied 
to crossing such higher-order resonances. But even if it 
can be applied, it is clear that the resonance strength 
cannot be obtained from a Fourier coefficient of uj in 
(|10|l . Moreover, in high energy accelerators, the n-th 
order Fourier coefficients of uj are not even the dominant 
contribution to the strength of a n-th order resonance. 
Usually the former contain G7, whereas the feed-up con- 
tributions from combining m lower order harmonics, con- 
tain (G"f) m , which can be an exceedingly large number. 

Only for first-order resonances, where J2t=i life I = 1> 
there is no feed-up contribution and the Fourier compo- 
nents can generally be used in the Froissart-Stora formula 
and there are different straight forward ways of comput- 
ing e K in that case 0,0- 



B. Solutions for the SRM 

The analytically solvable model advertised above is 
usually called the single resonance model (SRM). It has 
^0 = ^0^0 and an uj which only has one Fourier contribu- 
tion, uj = e tc (mcos$ + ^sin<i>), with $ = j + j ■ <f> + $ - 
Note that the modulus of its higher-order Fourier coeffi- 
cient is denoted as e K since there are no lower order coeffi- 
cients that could contribute to the resonance strength by 
a feed-up precess. Any dependence on the orbital actions 
can be expressed by e K (J). 

This uj is perpendicular to no and tilts spins away from 
Hq . Since 4* $ = Q, the frequency is k = jo + j ■ Q and 



the equation of motion (|10|l becomes 



_L S = ivos - ie K e^ e +^ V^W • (12) 

When the coordinates in the [m, i, no] system are ar- 
ranged in column vectors [Til ITa ] , one obtains 

1 j ( e K cos $ 

Initial coordinates z*i are taken into final coordinates Zf 
after one turn according to the relation <&/ = $^ + 2irQ 
whence $/ = $^ + 2itk. Now the orthogonal matrix 
T(e, <p) is introduced to describe a rotation around a unit 
vector e by an angle if. Transforming the spin com- 
ponents of S into a rotating frame using the relation 
Sr = T(n*o, — $) • 5, one obtains the simplified equation 
of spin motion 



— Sr — flu X Sr , Qr — 

du 




5 = v 



(14) 



If a spin field is oriented parallel to VLr in this frame, 
it does not change from turn to turn. Therefore ur — 
Qr/\Qr\ is an invariant spin field. In the original frame, 
this n-axis is 



1 



e K cos $ 



n($) = sig((5)- j e K sin$ 



(15) 



where the 'sign factor' sig(5) has been chosen so that 
on the closed orbit (e K = 0) the n-axis n(<I>) co- 
incides with no = (0, 0, 1) T . As required, n is 
both a solution of the T-BMT equation {T^Jl, f^ft = 

sig(<5) ■ i ^ is - ( — sin<!>, cos*! 1 , 0) T = O x n and, as with any 
function of phase space, a 27r-periodic function of the 
angle variables $ and of 9. 

This analytically solvable model can also be used to 
illustrate the construction of a phase independent but 
amplitude-dependent spin tune v{ J)- Once an n-axis 
has been obtained, one can transform the components 
of S into a coordinate system [ui,U2,n\. With the sim- 
ple choice 



w 2 ($) 

Ml($) 



no x n 



n x n 



= sig(<5) 




1 



S cos $ 



— I <5sm$ 
■ e K 



(16) 



(17) 



Mi is equal to 112 x n and the basis vectors are clearly 
27r-periodic in $ and in as required. Since ft and the 
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basis vectors Mi and U% comprise an orthogonal coordi- 
nate system for all 9, and since n precesses around f2, 
one has -§gu 2 = (fi — vri) x u 2 with the rotation rate v 
which can be computed by the relation 



( - n X U 2 ) • Ui 

— K COS $ + I^o cos $ 

sig(<5)[| — K sin $ + vq sin$ 



sig(<5)A 



(18) 



In general, the so found rotation could depend on <I> 
and an additional rotation of u\ and U2 around ft can now 
be used to make v independent of the angle variables $ 
and to define the amplitude-dependent spin tune. Here 
however, v is already independent of $ and it is therefore 
an amplitude dependent spin tune, and e K = |c3(z)| char- 
acterizes the orbital amplitude. The freedom of rotating 
Mi and U2 around n for each phase space point can be 
used to obtain a v which reduces to vq on the closed orbit 
(e K = 0). Wc let Mi and u 2 rotate around n by — $, to 
give the amplitude-dependent spin tune 



v — sig(£)A + k 



(19) 



The corresponding uniformly rotating basis vectors U\ 
and u 2 become 

u\ = ui cos$ — {t2sin<£> , u 2 — u 2 cos & + u\ sin $ . (20) 

On the closed orbit, the coordinate system now reduces 
to 

n — > n , mi — > sig(S)m , u 2 — > sig(<5)Z , v — > v Q . (21) 

This model leads to the average polarization on the torus 
with e K (J), 



mm 



A = V — K , 5 = l>{) — K , 



1 



f) , (22) 
(23) 



where the distance of the amplitude-dependent spin tune 
v from the resonance has been denoted by A, which is 
equivalent to sig((5)A. In Fig. (top) P\i m is plotted 
versus v$. It drops to at vq = k since according to 
(|T5|l the cone of vectors {n(<I>)|$ G [0, 2n]} opens up for 
small values of \S\. This strong reduction of P\i m occurs 
when v approaches n, i.e. close to spin-orbit resonances. 
According to (|19fl v is never exactly equal to k, but it 
jumps by 2e K across the resonance condition v = k, which 
is shown in Fig. ^ (bottom). This jump of the spin tune 
could in principle be transformed away since the sign of 
the spin tune depends on the sign of the rotation direction 
n. Here the sign of n in l|15(l has been fixed by choosing 
ri ■ n > 0, and the tune jump is therefore essential. 

Now we want to investigate the crossing of resonances, 
for the SRM, and describe spin motion when a parameter 
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FIG. 1: Pum and the amplitude-dependent spin tune v(e K ) for 
the SRM in the vicinity of vq = k, for k — 0.5 and e K =0.1 



r of the system is being slowly changed, i.e. -^t = a. In 
particular this allows the study of an acceleration where 
vq crosses the frequency k. It is useful to describe the 
spin motion in the coordinate system [mi, u 2 ,n]. In order 
to take account of the change of the basis vectors with 
the parameter r, we use that for a vector with \ui\ = 1, 
d T Ui is perpendicular to ui so that it can be written as a 
rotation, 



d 



d 



—n = rjxn , —ui 

OT OT 



rj x Mi . 



The rotation vector if is then given by, 

fj = — (mi x d T u 1 + u 2 x d T u 2 + n X d T n) 



(24) 



(25) 



Since v$ = G7 in a flat ring, the acceleration process 
in the SRM is usually described by a slowly changing 
i>a = k + t with t = a6 while assuming that k and e K 
do not change with energy. This leads to the following 
expressions for the variation of the basis vectors and for 
77: 



(26) 
(27) 
(28) 



d T u\ 


= sig(<5)-^r^cos$ 


d T u 2 


= sig(<5)-^r^sin$ 


d T n 


= -sig(£)-p«i , 


V 


= 8ig(*)^-(-U2- 







(29) 



In a general system, the equations of motion for the 
components of S = uisi + u 2 s 2 + nJs are described as 



ol{ 773S2 - mJs )-v{J,t)s 2 
a( ViJs ~ V3S1 ) + v(J,t)s 1 
a( T] 2 Si ~ r)is 2 ) 









■( 







(30) 
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In complex notation with s — Si + is2, rj = Tji + ir/2, and 
Js = \/l ~~ l^l 2 ; this gives 

-^-s = i[v{J,r) - ar] 3 ]s + iarjy/l - |s| 2 . (31) 

For the SRM, the Eqs. JS) and QJJ lead to 77 = 
_ i e te i(«e+*o) ; ^ = o, and 

±3 = z[sig(<5)A + k]s + a^^+^y/l^W ■ (32) 

Note again that the spin tune sig((5)A + K in this equation 
jumps by 2e K at = k. 

We will now describe how this equation for the SRM 
leads to the Froissart-Stora formula. After that, we will 
use the similarity of the SRM in Eq. i|32|) and the equation 
for a general system in Eq. (|31(l to show how higher-order 
resonance strength can be introduced and how they can 
be computed. 

III. THE FROISSART-STORA FORMULA 

For the analytically solvable SRM the change of the 
adiabatic invariant J5 — S ■ ft can be computed explic- 
itly. When the design-orbit spin tune changes during 
the acceleration process, resonances will be encountered, 
where v jumps from «±e K to re=Fe K while the spin is under 
the strong influence of an approximately resonant Fourier 
contribution of to. It is then found that for some speeds 
of the spin tune change, parametrized by a, a reduction 
of polarization can occur is due to a generally irreversible 
reduction of Js rather than a temporary decrease of Pum, 
and which does not recover after the energy has increased 
and the resonance is crossed. 

To describe the reduction of polarization during reso- 
nance crossing, (|32[) can be used but the usual approach 
is to insert a changing closed orbit spin tune vq into the 
equation of motion 1|13[) . The method of solution de- 
pends on the form of the function v (8) 0, GJ, El EI . If 
the closed-orbit spin tune changes like vq = n + a8, the 
equation of spin motion can be solved in terms of con- 
fluent hypergeometric functions. The equations for arbi- 
trary initial conditions are quite complicated but when 
at 9 — » —00 a vertical spin S3(— 00) = 1 is chosen as the 
initial condition then the vertical component at 8 — > +00 
is given by the well known and regularly used Froissart- 
Stora formula (l6| . 

s 3 (oo) = 2c~^^ - 1 . (33) 

In the case of a strong perturbation e K , or when the 
acceleration is very slow, spins follow the change of 
ri(Q). The 71-axis in l|15l) has a discontinuity from n_ = 
— e K (cos$,sin$,0) T just below resonance to ri + = — n_ 
just above resonance. Spins do not follow this instanta- 
neous change of sign, but they then follow — n adiabat- 
ically after the resonance has been crossed. Therefore 



S3 (00) is close to —1 for a slow change of uq. When the 
perturbation is weak or crossed very quickly, then spin 
motion is hardly affected and 53(00) is close to 1 in ij^l . 
In intermediate cases, | S3 1 is reduced. In the first case 
the polarization is preserved but the spins are reversed. 
In the second case the polarization is preserved without 
reversal. In the third case the polarization is no longer 
vertical but precesses around the vertical so that the time 
averaged polarization is reduced. 



IV. THE FROISSART-STORA FORMULA FOR 
HIGHER-ORDER RESONANCES 

As mentioned above, the Froissart-Stora formula in 
Eq. is regularly used to describe the reduction of 
polarization due to vertical betatron motion during reso- 
nance crossing in accelerators where the closed-orbit spin 
tune vo changes with energy. These descriptions were 
normally restricted to flat rings and v§ = G7. 

Since Siberian Snakes [HEl El EH are unavoid- 
able for high-energy polarized beam acceleration, the 
design-orbit spin tune is i in most cases which will be 
considered here and it does not change during accelera- 
tion. Since the orbital tunes are never chosen to be |, 
first-order resonances with v = jo ± Qk are avoided and 
higher-order resonances can become dominant. But since 
the strength of such resonances cannot be obtained as a 
Fourier coefficient of ui(z(6),6), a method for obtaining 
the strength of the higher-order resonances is required in 
order to use the Froissart-Stora formula when Siberian 
Snakes are in use. 

HERA-p will require at least 4 Siberian Snakes 0, H3, 
I2H Efll ] . The snake angles tpj of these 4 snakes can be 
chosen quite arbitrarily, except for the restriction Aip = 

— <P3 + f2 — <Pi = § • To illustrate crossing higher-order 
resonances a snake scheme for HERA-p was chosen that 
has 4 Siberian Snakes with snake angles of ?, 0, ? and 
in the South, East, North and West straight section, 
respectively. 

In Fig. [2] the amplitude-dependent spin tune (green) 
and P\i m (blue) are plotted versus the reference momen- 
tum for a vertical amplitude of IQtt mm mrad. Many 
higher-order resonances can be observed. The curves for 
P\i m and v{J) were computed with the non-perturbative 
algorithm SODOM II 1311 using the spin-orbit dynam- 
ics program SPRINT |l5l |32( . The n-axis and also P\i m 
are in general different at different azimuth 80 . For this 
figure and for all following plots of Pu m , the n-axis was 
observed at the interaction point of the ZEUS experiment 
in the South of HERA. 

While the design-orbit spin tune remains at |, the 
amplitude-dependent spin tune v{J y ) changes with en- 
ergy and is in resonance with 2Q y at the second line (red) 
and with 5Q y — 1 at the bottom line at several energies. 
In both cases a clear change of Pu m can be observed. The 
reduction of Pu m at some resonances is similar to the be- 
havior for the single resonance approximation shown in 
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where Pu m is reduced at those resonances. The drop 
of P Vlm at 811.2 GeV/c is due to the 2 — 5Q y resonance, 
which lies a little below the 2Q y line. At all other en- 
ergies where this resonance is crossed, no influence on 
P\im can be observed since the corresponding fifth-order 
resonance strength is very small. At some second-order 
resonances, Pum increases resonantly. Presumably, two 
resonant effects are in constructive interference at these 
energies. Nonetheless, polarization can be reduced when 
these resonance positions are crossed during acceleration 
since a sudden increase of P\i m ={n) is due to a sudden 
change of n(z) which might be too sudden for the adia- 
batic invariance of Js = S ■ n(z) to be maintained. In 
addition one can see in Fig. [3 that the spin tune v(J y ) 
has discontinuities at some of the resonances. 

When spin motion in a ring is approximated by a sin- 
gle resonance with k = jo ± Q y and then Siberian Snakes 
are included in the ring, it has often been noted that 
only odd-order resonances with k — jo + j y Q y appear, 
i.e. j y is odd. However it can be shown by nonlinear nor- 
mal form theory that this is a feature of any ring with 
midplanc symmetric spin-orbit motion and is not pecu- 
liar to rings with Siberian Snakes 0j . For rings without 
midplanc symmetry, resonances of even order can appear 
also. HERA-p has non-flat regions, and rings with closed- 
orbit distortions in general do not have midplane sym- 
metric motion. Then, resonances with even j y can also 
appear and be destructive. In fact, the resonances with 
j y = 2 are among the most destructive spin-orbit res- 
onances in HERA-p after Siberian Snakes are included. 
For the IUCF cooler ring with a partial snake running, 
second-order resonances have been observed experimen- 
tally m. 
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FIG. 2: Pum (blue) and v(J y ) (green) for particles with a 
4.2cr vertical amplitude of 70-7T mm mrad in HERA-p with 
and Q y — 0.289 . Three resonance lines cross v and at each 
crossing Pu m exhibits a large variation and there are jumps 
in v, bottom: v — 5Q y — 1, middle: v — 2 — 5Q y , and top 

v = 2Qy 

When a parameter r is being varied, the spin motion is 
described in the coordinate system [u\,U2,n] by Eq. (|31|1 . 
In the following we will demonstrate that this equation 
has some characteristics of the equation of spin motion 
(|32|l of the SRM. If the spin tune v has a discontinuity 



from v_ to v + at some energy, then we define the center 
frequency k* — \(v- +v + ). To take the jump of v into 
account, we introduce A* = \v — k*\, which does not 
have a discontinuity and we express the spin tune as v — 

SVg(v - K*)A* + K*. 

Since if is related to the basis vectors by (|25(l . it is a 
27r-periodic function of $ and 8. The jump of v across k* 
can be produced by a Fourier component of rj if there is a 
set of integers so that j ' • Q + jo = K*- This is the case in 
all instances of spin tune jumps presented here. Accord- 
ingly, one can analyze what happens when the Fourier 
component rj K *e^ K e +*o) of 77 dominates the motion of 
s. For that analysis, all other Fourier components of r\ 
are ignored. When a is small, spins which are initially 
almost parallel to the n-axis remain close to n so that s 
is small and ar) 3 s can therefore be ignored. This leads to 

= i(sig(p - k*)A* + k*)s + /ctt7 K .e iK * e+ *° \s\ 2 . 

w 

(34) 

Due to its similarity with (|32|l , this equation will produce 
the observed spin tune jump by 2e K * = \v + — v_ \ if i] K * = 
= (t/-K**) 2 m * ne v i cm ity of the energy where the 
jump occurs. Otherwise (|34J) would not reproduce this 
jump. One is then left with a relation which has exactly 
the structure of the equation of motion loL'l) for the SRM. 
Therefore, the Froissart-Stora formula can be applied to 
estimate how much polarization is lost when a polarized 
beam is accelerated through the energy region where the 
spin tune jumps by 2e K . In the following we will check 
whether, for some higher-order resonances in HERA-p, 
all assumptions leading to the approximation (|34J) are 
satisfied to the extent that the Froissart-Stora formula 
describes the reduction of polarization well. 

The basis vectors ft, Ux, and U2, and the amplitude- 
dependent spin tune v can in general only be computed 
by computationally intensive methods. The perturbing 
function r\ is then obtained from 

r\ = fj ■ (u\ + \U2) = ff ■ (—ft x H2 + in x ui) 

= (if x n) ■ (-u 2 + \u\) = i(ui + 1U2) • (d T n) ,(35) 

but the required differentiation is prone to numerical in- 
accuracies. However, when n is computed by perturba- 
tive normal form theory using differential algebra (DA) 
p34|. the differentiation with respect to r can be per- 
formed automatically. After rj is computed, the Fourier 
integral over the complete ring would finally be required 
in order to compute e K . 

If 131f) can be approximated well by a SRM, there is 
however a different and much less cumbersome method 
for determining the relevant resonance strength and 
the resonant frequency. Observation of the amplitude- 
dependent spin tune v( J) allows the determination of all 
parameters which are required to evaluate the Froissart- 
Stora formula for higher-order resonances: The spin tune 
jumps by 2e K , the center of the jump is located at the fre- 
quency k itself, and the rate of change of v with changing 
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energy is used to determine the parameter a for i|33|l . In 
the SRM this parameter is v °g K where i/ is the frequency 
of spin rotations when the resonance strength vanishes. 
Here the corresponding frequency, which would be ob- 
served if no perturbation 77 were present, is not directly 
computed. But it can be approximately inferred from 
the slope d T v at some distance from the resonance. 
According to (J2HJ, (n) is given by P^ M = 
/± — in the SRM. To check whether the observed 

drop of Pnm indeed shows the characteristics of the SRM, 
the width of the resonance dip in P^ M was obtained 
from the amplitude-dependent spin tune alone and then 
compared to the width of the dip in the actual Pu m of 
the system. This analysis was done for HERA-p's reso- 
nance at approximately 812.4 GeV/c and the results are 
shown in Fig. |3| The top left plot shows the dependence 
of P\im and v on the reference momentum for a vertical 
amplitude of 70n mm mrad which, with HERA-p's cur- 
rent one sigma emittance of 4tt mm mrad, corresponds 
to the amplitude of a 4.2<r vertical emittance. The mo- 
mentum range is as in Fig. [21 The low P\ im shows that 
many perturbing effects interfere in this region. In units 
of 7r mm mrad, the vertical amplitude of the particles 
in the top left graph is 70, in the middle graphs it is 40 
and 60, and in the bottom graphs 80 and 100. The hori- 
zontal scale displays the distance Ap in GeV/c from the 
momentum at the resonance. 

In the 4 bottom graphs, Pu m and P^ M are plotted 
for different orbital amplitudes, and the different reso- 
nance strengths are obtained from the jump in v{J y ). 
Only information about v was used to compute P^^ M ■ 
To allow better comparison, a linear change of P^ M 
with momentum was added as a background curve and 
the height of the dip was scaled to fit the actual Pu m . 
The width however was not changed. The distance be- 
tween spin tune and resonance has been magnified by 10, 
v* = k + 10(k — v) in these graphs. The tune jump is 
symmetric around the resonance line v — 2Q y , showing 
that a second-order resonance is excited. 

As shown in Fig. (top right) the tune jump scales 
approximately linearly with the orbital action variable 



with the actual drop in Pi 



J y . This is consistent with the crossing of a second-order 
resonance, since a frequency of 2Q y can be produced by 
monomials of ^JH y 'e ±ll ^ ve with order larger or equal to 
2. This linear scaling is not exact for two reasons: (1) 
The jump does not reduce to at J y = but already at 
some finite amplitude at which v(J y ) does not cross the 
resonance line. (2) When the amplitude is changed, the 
momentum at which the resonance occurs changes, and 
the resonance strength is in general different at different 
energies. Deviations from a linear dependence should 
therefore be expected. Pn m is already very low away 
from the resonance at v = 2Q yi indicating that other 
strong perturbations distort the invariant spin field and 
can interfere with the resonance harmonic. 

Thus we conclude that the resonance width computed 
in terms of the tune jump 2e K agrees surprisingly well 
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FIG. 3: Top left: Pum and v in the vicinity of the reso- 
nance at approximately 812.4 GeV/c for HERA-p. The dis- 
tance between v and resonance has been magnified by 10, 
v* = k + 10(k — v). Top right: Proportionality between 
tune jump 2e^ and the amplitude 2J y of a vertical emit- 
tance. Middle and bottom: Correlation between the width 
of the actual drop of Pn m and the predictions of the single 
resonance approximation using only the amplitude-dependent 
spin tune. Vertical amplitudes of particles in HERA-p in units 
of 7r mm mrad from top left to bottom right: 70, 40, 60. 80, 
and 100. Ap: distance from the momentum at resonance in 
GeV/c 

Since the higher-order resonances analyzed here show 
the established and characteristic relation between tune 
jump and reduction of Pi™, the applicability of the 
Froissart-Stora formula will now be tested. 

In Fig. 0| (top) Pij m and v are shown for HERA-p. Pn m 
is reduced at two resonances with v = 2Q y . The vertical 
tune had been chosen as Q y = 0.2725 so that these res- 
onances are crossed already for the small 0.75cr vertical 
amplitude of 2.257T mm mrad. At this small amplitude 
P\im is reasonably large. 

The spins of a set of particles were set parallel to the 
invariant spin field n(z) so that all had J5 = 1 at the 
momentum of 801 GeV/c. The n-axis had been com- 
puted by stroboscopic averaging Due to the rather 
large Pn m at that energy the initial polarization was ap- 
proximately 97%. Starting with this spin configuration, 
the beam was accelerated to 804 GeV/c at various rates. 
The average (Js)n over the tracked particles is plotted 
versus acceleration rate in Fig. 0] (bottom) together with 
the prediction of the Froissart-Stora formula. The av- 
erage (Js)n describes the degree of beam polarization 
which could be recovered due to the adiabatic invariance 
of Jg when moving into an energy regime where n{z) is 
close to parallel to the vertical. 
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The resonance strength €2Q y has been determined from 
the tune jump. The parameter a is proportional to the 
energy increase per turn cIe and is determined from the 
tune slope jg in Fig. 21 (top right) by the relation a = 

1 Av j 

The polarization obtained by accelerating particles 
through the second-order resonance agrees remarkably 
well with the Froissart-Stora formula. For the slow accel- 
eration of about 50 keV per turn in HERA-p, the polar- 
ization would be completely reversed on the 0.75 sigma 
invariant torus. This would lead to a net reduction of 
beam polarization, since the spins in the center of the 
beam are not reversed. 
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FIG. 4: Top: Pu m and v for a second-order resonance of 
HERA-p with Q y = 0.2725 and a 0.75<r vertical amplitude of 
2.257T mm mrad. Bottom: {Js)n after acceleration from 
801 GeV/c to 804 GeV/c with different acceleration rates 
(blue points) and the prediction of the Froissart-Stora for- 
mula (red curve) using parameters tiQ y and a obtained from 
u 

This result on the applicability of (|33() for the reso- 
nance strength and a obtained from the amplitude de- 
pendent spin tune is so important for detailed analysis of 
the acceleration process that it will be checked in another 



case. In the next example, the same lattice is used, the 
tune was adjusted to a realistic value of Q y = 0.289 and 
a 4.2cr vertical amplitude of 70-7T mm mrad was chosen. 
At this large amplitude, the second and fifth-order res- 
onances already shown in Fig. [21 are observed. Particles 
were then accelerated from 812.2 GeV/c to 812.6 GeV/c 
with different acceleration rates. Note that the initial 
condition has a vertical polarization of only 60%. Nev- 
ertheless this state of polarization corresponds to a com- 
pletely polarized beam, and 100% polarization can poten- 
tially be recovered by changing the energy adiabatically 
into a region where ft(z) is tightly bundled. These stud- 
ies emphasize again the importance of choosing n(z) as 
the initial spin direction. For example if the spins were 
initially polarized vertically, they would rotate around 
n(z) and that would lead to a fluctuating polarization, 
even without a resonance and it would not be possible to 
establish a Froissart-Stora formula for higher-order reso- 
nances. 

As shown in Fig. EI P\im is as l° w as 0.11 in the center 
of the displayed region. Obviously other strong effects 
beyond the second-order resonance are present and over- 
lap with it. The bottom figure shows (Js)n after the 
acceleration. The fact that {Js)n is again described very 
well by the Froissart-Stora formula ll-i-il) is an impressive 
confirmation of the conjecture. 

The two data points at largest acceleration speed in 
Fig. 0] (bottom) are lower than predicted by the Froissart- 
Stora formula. A possible explanation is the following: 
at very large acceleration speeds the resonance region is 
crossed so quickly that the spin motion is hardly dis- 
turbed. But when the n-axis n_ before the resonance 
region is not parallel to the n-axis ri + after the resonance 
region, then the spins which initially had Jg = 1 will ap- 
proximately have Js — n- • n + after the resonance region 
is crossed, which is smaller than the Froissart-Stora pre- 
diction, which approaches 1 for large acceleration speeds. 

Here the parameter t was the slowly changing momen- 
tum. This generalized way of using the Froissart-Stora 
formula can however also be used when other system pa- 
rameters change. An example can be found in |35J, where 
the particle's phase space amplitude is changed artifi- 
cially slowly in order to compute the invariant spin field 
at various orbital amplitudes. In 10] an example is dis- 
played where the Froissart-Stora formula is successfully 
applied to a resonance which is encountered because of a 
slow variation of Q y . 



V. THE CHOICE OF ORBITAL TUNES 

When the amplitude-dependent spin tune v(J) of par- 
ticles with the amplitude J crosses a resonance, for exam- 
ple during acceleration, the beam polarization is usually 
reduced. It is therefore important to find suitable orbital 
tunes so that low-order spin-orbit resonances are far away 
from the operating point. In particular, when Siberian 
Snakes are used to maintain a closed orbit spin tune of \ , 
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FIG. 5: Top: Pum and v for a second-order resonance of 
HERA-p with realistic tune of Q y = 0.289 and a large 4.2a 
vertical amplitude of 70n mm mrad. Bottom: {Js)n after 
acceleration from 812.2 GeV/c to 812.6 GeV/c with differ- 
ent acceleration rates (blue points) and the prediction of the 
Froissart-Stora formula (red curve) using parameters eiQ y and 
a obtained from v 



it is important that these Snakes are optimized so that 
higher-order resonances do not lead to large deviations 
of the amplitude dependent spin tune from this value. 
Such optimal choices of snakes are discussed in [2|j . The 
dominant effects are due to radial fields on vertical beta- 
tron trajectories. Thus Fig. (right) shows the resonance 
lines v = jo + jQ y up to order 10 in the v-Q y plane. If 
the spin tune on the closed orbit is fixed to Uq = | by 
Siberian Snakes the orbital tune can be chosen to avoid 



resonance lines. However, the dynamic aperture of pro- 
ton motion should not be reduced and the tunes have to 
be far away from low order orbital resonances. Figure 
(|rj|) (left) shows the Q x -Q y tune diagram with resonance 
lines up to order 11. The operating point has to stay 
away from these resonance lines. 

The established tunes of HERA-p operation Q x = 
0.294, Q y = 0.298 or Q x = 0.298, Q y = 0.294 (red 
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FIG. 6: Left: the current orbit tunes (Q x = 0.294,Q H = 
0.298) or (Q x = 0.298,Q H = 0.294) (red) and the new orbit 
tunes for polarized proton operation (Q x — 0.291,Q H = 0.289) 
or (Q x = 0.2675, Q y = 0.271) (blue) in the x-y resonance 
diagram. All resonances up to order 11 are shown. Differ- 
ence resonances are indicated in green. Right: The current 
vertical tunes (red) and the new vertical tunes (blue) in the 
spin-orbit resonance diagram. The odd spin-orbit resonances 
(black) and the even spin-orbit resonances (green) are shown 
up to order 10 in the vicinity of closed-orbit spin tune i/q = i. 
For HERA-p, the resonances of second order (fat green) and 
of fifth order (fat black) are most destructive 



points) would be unfortunate choices due to their close- 
ness to the resonance v = jo ± 5Q y . For HERA-p with 
Siberian Snakes, several simulations have shown that the 
resonances of second order and of fifth order are most 
destructive. This is supported by Fig. |3 Therefore two 
new tunes (blue points) are suggested which have an opti- 
mal distance from low-order spin-orbit resonances. It has 
been tested experimentally that HERA-p could operate 
at these tunes. 
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